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SUMMARY 


Simultaneous  optical  and  meteorological  measurements  between  ground  level  and  6  km  altitude  gathered  by  the 
Visibility  Laboratory  in  recent  years  in  the  United  States  and  Europe  with  an  instrumented  aircraft  have  provided  the  basis 
for  the  development  of  a  technique  for  estimating  atmospheric  path  radiance  and  directional  contrast  transmittance.  The 
technique  is  relatively  fast  and  easy  to  apply  yet  is  designed  to  take  full  advantage  of  the  present  limited  capability  to 
observe  and  predict  the  relevant  atmospheric  variables.  Various  modelling  approximations  are  introduced  to  simplify  the 
calculations  of  radiative  transfer  and  to  help  specify  the  vertical  distribution  of  atmospheric  attenuation  coefficient  and  the 
single  scattering  phase  function.  Critical  examinations  of  the  assumptions  and  approximations  were  made  through 
reference  to  the  extensive  series  of  airborne  meteorological/ optical  measurements.  Model  input  parameters  are  the 
number  of  atmospheric  layers  selected  and  the  average  optical  scattering  ratio  and  single  scattering  albedo  for  each  layer, 
the  solar  zenith  angle,  the  extra-terrestrial  solar  irradiance.  a  representative  wavelength  and  the  surface  reflectance.  The 
computer  program  calculates  the  sky  and  terrain  radiance  corresponding  to  any  selected  point  and  direction  and  the  contrast 
transmittance  for  any  slant  atmospheric  path. 

The  changes  in  image  transmittance  characteristics  that  are  associated  with  the  natural  variations  in  environmental 
physical  parameters  are  examined  through  systematic  application  of  the  modelling  techniques.  Comparisons  are  made  of 
the  calculated  changes  in  visible  spectrum  contrast  transmittance  that  result  from  typical  changes  in  the  optical  thickness  of 
the  low-level  haze  layer,  the  aerosol  absorption,  the  surface  spectral  reflectance,  the  solar  zenith  angle  and  the  viewing 
path 
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AN  OPERATIONAL  TECHNIQUE  FOR  ESTIMATING 
VISIBLE  SPECTRUM  CONTRAST  TRANSMITTANCE 


1.  INTRODUCTION 

The  visual  detection  of  distant  objects  is  restricted 
by  air  molecule  and  aerosol  particle  scattering  and  to  a 
lesser  extent  by  atmospheric  absorption.  The  complexity 
of  multiple  scattering  and  absorption  processes  and  the 
inherent  variability  of  the  aerosol  particle  distribution  and 
composition  require  effective  simplifications  in  order  that 
estimates  of  contrast  transmittance  through  the  atmo¬ 
sphere  can  be  made  rapidly  and  consistently. 

A  series  of  modelling  approximations  that  relate 
optical  scattering  properties  to  meteorological  variables  are 
derived  from  experimental  data  and  combined  with  avail¬ 
able  analytic  approximations  for  radiative  transfer  calcula¬ 
tions  to  develop  an  operational  technique  for  the  estima¬ 
tion  of  directional  path  radiance  and  directional  slant-path 
contrast  transmittance.  The  approach  involves  validation 
of  the  individual  elements  of  the  procedure  through  refer¬ 
ence  to  simultaneous  optical  and  meteorological  measure¬ 
ments  gathered  by  the  Visibility  Laboratory  (Johnson 
etal..  1979)  over  the  past  several  years  with  a  specially 
instrumented  Air  Force  C-130  aircraft  under  sponsorship 
of  the  Air  Force  Geophysics  Laboratory  The  high  resolu¬ 
tion  profile  measurements  were  made  in  a  wide  range  of 
meteorological  and  geographical  conditions  in  the 
United  States  and  western  Europe. 

The  path  contrast  transmittance  in  the  visible  spec- 
<rum  depends  most  critically  upon  the  distribution  of  total 
volume  scattering  coefficient,  both  along  the  path  and  in 
the  surrounding  atmosphere.  Accordingly,  an  attempt  was 
made  <o  condition  the  approximation  procedures  for  the 
radiative  transfer  calculations  and  for  all  input  parameters 
in  a  manner  consistent  with  the  existing  capability  to 
model  and  predict  the  scattering  coefficient  structure  and 
behavior 


2.  THEORETICAL  DEVELOPMENT  AND 
COMPARISONS  WITH  OTHER  MODEL  RESULTS 

Neglecting  turbulence  effects,  the  equation  for  the 
apparent  spectral  radiance  of  the  background  b  at  range  r 
along  the  path  of  sight  specified  by  zenith  angle  8  and 
azimuthal  angle  <b.  can  be  written  (Duntley  ei  al..  1957) 

hL,(:,fi.<b)  -  r,(z,»)bLJ:,.8.<t>)  +  Z.r'<z.».*>  .  (I) 


where  bL0  is  the  inherent  background  radiance  at  target 
altitude  z„  T,  -  exp  -J  a(r)dr  is  the  path  transmittance,  a 

O 

is  the  volume  attenuation  coefficient,  and  L'  is  the  path 
radiance  produced  by  the  scattering  of  light  from  the  sun 
and  from  the  surrounding  sky  and  terrain  into  the  path  of 
sight. 

The  path  radiance  is  given  by 

r 

Oz.8.6)  -  f  LA:‘.9,<b)  T,h'.8)dr'  .  (2) 


where  LJ:'.9.d>)  is  the  path  function,  defined  as  the  point 
function  component  of  path  radiance  generated  by  the 
directional  scattering  of  light  reaching  that  point  of  the 
path.  The  expression  for  the  path  function  can  be  written 
in  terms  of  its  contributions  from  the  scattering  of  solar 
scalar  irradiance  ,«(;)  and  from  the  scattering  of  sky  and 
earth  radiances  L(z. »'.<*'•  as  follows  (Gordon,  1969): 

L.l:.8.<b)  -,<(;)<  r<:./3,)  +  J  Ui.8'.<b‘)<r(z,p)dll' ,  (3) 

4»r 

where  <r(z,0)  is  the  directional  volume  scattering  function 
at  angle  /S  between  the  path  of  sight  and  direction  of  the 
source  light,  d n  is  an  element  of  solid  angle,  and  lis  is  the 
scattering  angle  with  respect  to  the  sun. 

In  turn,  the  spectral  contrast  transmittance  of  the 
path  of  sight  can  be  expressed  (see  Duntley  etal.,  1957) 
directly  as  the  product  of  the  path  transmittance.  Tr(z.fl), 
and  the  ratio  of  the  inherent,  bL„,  and  apparent,  bL,, 
background  radiances  as  follows: 

Ck.*.i)/f„(z„W)-  Tr{i.9)iLc(z,.«.*)/tLr(:.8.4>)  .  (4) 

where  C,  -  (,L,  -  bLr)/bL,  is  the  apparent  target  contrast 
at  path  length  r ,  C„  -  (,T„  -  bL„)/bL„  is  the  inherent  tar¬ 
get  contrast  at  altitude  z,  and  ,L„  is  the  inherent  target 
radiance.  Thus,  the  contrast  transmittance  of  a  given  path 
does  not  depend  upon  intrinsic  target  characteristics  but  is 
a  function  only  of  the  directional  radiance  distribution  in 
the  atmosphere  and  the  path  transmittance.  The  expres- 


sion  is  strictly  applicable  only  for  monochromatic  radiation 
but  may  be  applied  with  good  approximation  to  reasonably 
broad  spectral  bands  in  the  visible  portion  of  the  spec¬ 
trum. 

From  Eq.  (3)  we  note  that  the  path  function  and  in 
turn  the  contrast  transmittance  as  calculated  from  Eq.  (4) 
depend  upon  the  direction  of  the  viewing  path  relative  to 
the  distribution  of  light  reaching  the  path.  The  first  term 
on  the  right  hand  side  of  Eq.  (3)  is  the  contribution  of  pri¬ 
mary  scattering  of  direct  solar  irradiance.  The  directional 
volume  scattering  function  may  be  expressed 

<r(z,/3)  -  P(z,f})s(z)  ,  (5) 

where  s(z)  is  the  total  volume  scattering  coefficient,  and 
P(z,$)  is  the  single  scattering  phase  function,  which 
defines  the  probability  that  incident  radiation  will  be  scat¬ 
tered  in  the  direction  given  by  scattering  angle  0.  The 
phase  function  varies  significantly  with  the  scattering  pro¬ 
perties  of  the  aerosol  particle  distribution  in  the  atmo¬ 
sphere.  Modelling  procedures  used  herein  either  accept 
P(z,i 8)  as  calculated  from  theory  from  an  assumed  or 
measured  particle  size  distribution  or  use  an  analytic 
representation  of  P(z.p)  as  described  in  Section  4  of  this 
paper. 

Substituting  Eq.  (5)  in  Eq.  (3),  the  expression  for 
L.(z, »,</>)  becomes 

L.(z,9,<i>)  -  st(z)P(z,0)z(z)  (6) 

+  ff  L(z,8',<t>')P(z,p)s(z)dCl‘  . 


for  fast  calculation  of  the  path  radiance  component  due  to 
scattering  of  the  background  sky  and  terrain  radiances  at 
any  point  and  direction  provided  that  the  asymmetric 
influence  of  the  prominent  forward  scatter  peak  and  the 
irradiance  profile  are  managed  adequately.  The  delta- 
Eddington  approximation  introduced  by  Joseph,  Wiscombe 
and  Weinman  (1976)  satisfies  the  requirement.  It  differs 
from  the  standard  Eddington  approximation,  which 
assumes  a  simple  cosine  dependence  of  the  single  scatter¬ 
ing  phase  function,  in  that  it  approximates  the  phase  func¬ 
tion  by  a  truncated  forward  scatter  peak  and  a  two-term 
phase  function  expansion, 

4ir  Pd(fl)  -  If  8U  -  cos/3)  +  (1  -  /'Hi  +  3g'cos/3)  (7) 

where  /’  is  the  fractional  scattering  represented  by  the  for¬ 
ward  peak  and  g'  is  the  asymmetry  factor  of  the  truncated 
phase  function.  In  effect,  the  delta- Eddington  approxima¬ 
tion  transforms  most  of  the  enhanced  radiance  in  the  solar 
aureole  into  the  direct  solar  flux  component,  and  assumes 

/'(.’)  -  gHz)  .  (8) 

Joseph  et  al..  (1976)  show  that  calculations  of  radiative 
transfer  with  the  detta-Eddington  approximation  can  be 
carried  out  with  the  standard  Eddington  computer  code 
(Shettle  and  Weinman,  1970)  with  the  following  changes 
of  variable 

Ar'  *>  (1  -  u>f  )  At  .  (9) 


<*>’  ■>  0>(1  —  /')/(l  —  U»/')  ,  (10) 

Modelling  techniques  for  the  calculation  of  the  diffuse 
component  of  the  path  function  are  discussed  in  the  fol-  and 

lowing  paragraphs.  g,  _  {^nnx_n ,  (11) 


2.1  Calculation  of  the  Path  Function 

The  second  term  on  the  right  hand  side  of  Eq.  (6)  is 
the  component  of  the  path  function  resulting  from  the 
scattering  of  diffuse  radiance  reaching  the  path  from  the 
surrounding  sky  and  terrain.  It  has  a  directional  depen¬ 
dence  due  to  the  asymmetry  in  the  sky  and  earth  radiances 
illuminating  the  path  and  in  the  scattering  phase  function. 
Precise  numerical  calculation  of  the  path  radiance  resulting 
from  the  complex  multiple  scattering  processes  requires 
large  amounts  of  computer  time.  For  this  reason,  rapid 
approximate  methods  are  employed  extensively  for  radia¬ 
tion  transfer  calculations.  The  appropriate  choice  of  com¬ 
putational  method  from  among  the  variety  of  available 
methods  depends  upon  the  results  desired  for  the  applica¬ 
tion  at  hand. 

While  it  is  important  to  retain  complete  directional¬ 
ity  for  calculation  of  the  path  radiance  component  due  to 
single  scattering  of  direct  solar  radiance,  approximate 
hemispherical  two  stream  methods  can  be  used  effectively 


where  At  is  the  optical  thickness  of  the  layer  and  w  is  the 
single  scattering  albedo. 

As  an  integral  part  of  the  technique  for  estimating 
directional  contrast  transmittance,  the  products  of  the 
Eddington  computer  code  (Shettle  and  Weinman,  1970) 
are  used  directly  to  calculate  the  second  term  on  the  right 
hand  side  of  Eq.  (6). 

From  the  standard  Eddington  approximation,  the 
diffuse  radiance  is  assumed  to  be  given  by 

L(z,8‘,<b')  **  Z-£>(r)  +  Ld(z)cos0‘  .  (12) 


In  accordance  with  the  derivation  given  in  Appendix  A.  if 
we  substitute  Eq.  (12)  and  the  delta-Eddington  approxima¬ 
tions  given  by  Eq.  (7).  (8)  and  (11)  into  the  last  term  of 
Eq.  (6)  and  integrate  over  O'  and  <t>'.  we  have 
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//  z.(z,#',*)/>(.’./su(z)</n  (13) 

4?r 

-  s(zl  UD(z)  +  glD(zkos0l  . 

Commensurate  with  the  delta-Eddington  approximations 
and  Eq.  (9),  the  expression  for  the  solar  scalar  irradiance 
at  altitude  :  is 

,«(;)  -  ,«(!»>  expt-r'/cos#,)  .  (14) 

where  ,<(<»>  is  the  extraterrestrial  solar  scalar  irradiance. 

Substituting  Eqs.  (13)  and  (14)  into  Eq.  (6).  the  expres¬ 
sion  for  the  path  function  becomes 

LAz.Q.<t>)  -  s(;l  I  )  1«(“-)exp(-r'/cos«1 )  (15) 

+  L0  +  gLD< os«  I 

2.2  Calculation  of  Sky  and  Terrain  Radiance 

For  an  assumed  plane  parallel  and  horizontally 
homogeneous  atmosphere,  the  directional  path  radiance. 
L'(z.$.<b ),  and  the  radiance,  Hz. 0.6),  can  now  be  calcu¬ 
lated  from  Eqs.  (1),  (2)  and  (15)  through  finite  summa¬ 
tion  over  adjacent  atmospheric  layers  using  the  trapezoidal 
rule.  For  upward  paths  of  sight  the  inherent  background 
radiance  at  target  altitude  is  given  by 

bL0(z, .0,6)  ~  L’x  (z,,0,6)  ,  (16) 

where  Ll,  ( z,.0.6 )  is  the  path  radiance  as  determined  at  alti¬ 
tude  z,  for  the  path  of  sight  from  the  top  of  the  atmo¬ 
sphere  to  z,.  Looking  downward,  the  inherent  background 
radiance  is  given  by 

„LJz,.0.6)  -  Tm)z,,0)  R{H,6)E(Q.d)l  n  (17) 

+  L‘„(z,.0.6)  . 

where  R)0.6)  is  the  directional  surface  reflectance. 
f(O.rf)  -  .  TtOkosfl,  f  irlLfltO)  +  2/3Lo  (0)  I  is  the  downwel- 
ling  irradiance  at  the  surface,  and  T„,  is  the  transmittance 
of  the  st  int  path  from  the  surface  to  target  altitude 

2.3  Comparisons  with  Other  Model  Calculations 

The  extensive  series  of  optical  and  meteorological 
'  icasu  remen  is.  referenced  in  Section  1.  has  been  used  to 
test  individual  modelling  approximations  and  to  help  vali¬ 
date  the  integrated  application  of  the  techniques  For 
example,  in  Section  5  -del  calculations  of  sky  and  ter¬ 
rain  radiance  based  upon  the  procedures  set  forth  n  Sec¬ 
tions  2,  3  and  4,  are  compared  with  the  high  resolution 
measurements  of  the  spectral  radiance  made  with  the 
instrumented  aircraft 


Additional  tests  were  carried  out  to  evaluate 
specifically  the  approximations  for  the  radiative  transfer 
calculations  as  described  above.  Model  calculations  of 
radiance  distributions  were  compared  with  the  theoretical 
results  as  calculated  from  comprehensive  radiative  transfer 
mode's,  which  are  mathematically  precise  but  in  general 
require  relatively  large  computer  capacity  and  time  for  the 
radi-'rce  calculations.  For  these  tests,  the  comparative 
accuracy  was  determined  under  the  assumption  that  the 
physical  system  is  perfectly  observed,  and  the  input  data 
are  identical  for  the  approximate  operational  model  and 
the  comprehensive  numerical  models.  The  set  of  exam¬ 
ples  selected  for  tests  by  the  ad  hoc  Working  Group  of  the 
Radiation  Commission,  International  Association  of 
Meteorology  and  Atmospheric  Physics  (Lenoble.  1977) 
were  used  for  the  analysis.  The  results  of  the  spherical 
harmonic  method  listed  in  the  Working  Group  report  were 
chosen  for  a  reference  standard  since  comparative  calcula¬ 
tions  using  this  method  were  presented  for  all  five  test 
cases.  A  single  layer,  plane  parallel  and  homogeneous 
atmosphere  with  a  black  underlying  surface  was  prescribed 
in  each  case.  Three  cases  have  a  Haze  L  (Diermendjian, 
1969)  aerosol  particle  distribution  with  an  optical  thickness 
of  1,  and  the  other  two  have  assumed  cloud  particle  distri¬ 
butions  with  an  optical  thickness  of  64  The  phase  func¬ 
tion  for  single  scattering  and  the  single  scattering  albedo 
are  specified  for  each  case.  The  comparative  results  of  the 
radiance  calculations  as  summarized  over  all  listed  interior 
verification  levels  (5)  and  zenith  view  angles  (11)  are 
given  for  each  case  (including  a  subset  for  each  of  3 
azimuthal  angles  for  case  3)  in  Table  2.1,  and  some  graph¬ 
ical  comparisons  of  test  results  are  shown  in  Figs  2-1,  2-2, 
and  2-3. 

It  is  important  to  note  that  the  agreement  among 
the  more  rigorous  numerical  methods,  including  the 
Monte  Carlo  techniques,  as  given  in  the  Working  Group 
report,  was  generally  within  a  few  percent  depending  upon 
the  completeness  of  scattering  history  retained  in  the  cal¬ 
culations.  We  see  from  Table  2.1  and  Fig.  2-1  that  the 
departures  of  the  simplified  model  radiance  calculations 
from  the  reference  values  for  the  haze  cases  are 
significantly  larger,  yet  the  indicated  accuracy  is  good  when 
considered  in  context  with  the  uncertainties  and  complete¬ 
ness  of  the  observational  data  base  from  which  the  esti¬ 
mates  must  be  derived  in  operational  practice.  The 
simplified  model  requires  only  a  few  seconds  of  large  scale 
computer  time  for  each  case,  and  the  model  is  designed 
for  easy  microprocessor  application  for  real  time  require¬ 
ments  Except  for  the  subset  of  case  3  with  the  azimuthal 
view  angle  directly  upsun.  <t>  *»  0°,  the  logarithmic  rms 
error  of  specification  for  the  haze  cases  was  between  O  il 
and  0  16  (or  approximately  II  to  16  percent  rms  error) 
The  approximations  inherent  in  the  operational  model 
serve  to  smooth  out  the  azimuthal  asymmetry  in  the  cal¬ 
culated  radiance  component  due  to  the  scattering  of 
incident  diffuse  light,  hence  the  increase  in  rms  error  for 
I  he  subset  <j>  =  0°  It  is  expected  that  the  error  of  esti¬ 
mate  would  continue  to  increase  in  the  solar  azimuthal 
plane  in  the  forward  scattering  direction  as  the  solar  zenith 
angle  increases  beyond  the  assumed  60°  for  case  3 
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l(r)  -  sKU)  +  sM(z)  -  u(z)a(z)  , 


(19) 


The  results  for  the  cloud  cases  listed  in  Table  2.1 
and  shown  in  Figs.  2-2  and  2-3  are  based  on  the  asymp¬ 
totic  solution  to  the  model  equations,  where  the  radiance 
distribution  at  any  point  interior  to  the  cloud  is  given  sim¬ 
ply  by  the  equilibrium  radiance,  L„(z,9,4)  -  L.(z,9,<t>)/a(z). 
Although  the  technique  was  not  developed  for  a  cloud 
environment,  it  was  found  that  the  calculated  equilibrium 
radiance  represents  the  directional  distributions  of  spectral 
radiance  for  the  two  reference  cloud  layers  with  an  over¬ 
head  sun  with  good  accuracy.  However,  equilibrium  radi¬ 
ance  does  not  obtain  for  the  viewing  paths  directly  upsun 
and  the  asymptotic  solution  yields  excessively  large  values 
of  calculated  radiance  for  6-6,  -  0.  The  predominant 
error  anomalies  for  the  directional  radiance  L(r,0,0)  are 
not  included  in  the  comparative  statistics  for  the  cloud 
cases,  so  that  the  summaries  listed  in  Table  2.1  for  cases  4 
(no  absorption)  and  3  (with  absorption)  cover  10  zenith 
viewing  angles  extending  from  37a  to  180s  and  5  optical 
thickness  levels  within  the  cloud.  The  calculated  radiance 
reduces  by  several  orders  of  magnitude  at  optical  depth  32 
for  the  absorption  case  (see  Fig.  2-3)  as  compared  with  a 
much  smaller  decrease  with  increasing  optical  depth  for 
the  conservative  case  (Fig.  2-2).  The  equilibrium  radiance 
solution  provides  reasonably  good  specification  of  absolute 
radiance  at  large  optical  depths  for  both  cases  but  the  per¬ 
centage  departure  from  the  reference  model  calculations 
increases  with  increasing  depth  for  the  cloud  case  with 
absorption  (Fig.  2-3). 


3.  SPECIFICATION  OF  THE  VERTICAL 
DISTRIBUTION  OF  TOTAL  VOLUME 
SCATTERING  COEFFICIENT 

The  spatial  distribution  of  total  volume  scattering 
coefficient  is  the  major  determinant  of  visible  spectrum 
contrast  transmittance.  Techniques  for  the  specification  of 
the  scattering  coefficient  profile  have  been  investigated 
during  the  course  of  the  aircraft  measurement  and  analysis 
program  (Johnson  elal.,  1979).  Some  of  the  results  of 
these  studies  and  their  application  to  operational  modelling 
procedures  are  summarized  in  this  section. 

3.1  Scattering  Ratio  Specification 

For  profile  modelling  purposes,  it  is  important  to 
consider  a  conservative  measure  of  scattering  coefficient 
that  in  the  absence  of  local  aerosol  particle  sources  or 
sinks  does  not  change  appreciably  following  the  air 
motion.  The  optical  scattering  mixing  ratio.  Q(z).  is  such 
a  parameter.  As  the  vertical  mixing  within  an  identifiable 
atmospheric  layer  becomes  more  complete  £Mr)  becomes 
more  constant  with  height  within  the  layer. 

The  optical  scattering  ratio  is  defined 

(Hz)  -  j(z)/j»<z)  .  (18) 

where  v*(z)  is  the  total  volume  coefficient  for  Rayleigh 
scattering  at  altitude  z.  Note  that 


where  a(z)  is  the  extinction  coefficient  and  jv(z)  is  the 
aerosol  scattering  coefficient,  also  referred  to  as  the  Mie 
scattering  coefficient.  It  follows  that  the  Mie  scattering 
ratio  is  given  by 

$M(r)/i#(r)  -  Q(z)  -  1  (20) 

The  Mie  scattering  ratio  also  would  be  constant  under 
conditions  of  complete  aerosol  mixing. 

Profiles  of  Q(z)  derived  from  the  extensive  series  of 
airborne  optical  measurements  made  by  the  Visibility 
Laboratory,  reveal  large  variability  depending  upon  the 
aerosol  particle  source  strength  and  the  nature  of  the  con¬ 
vective  and  turbulent  mixing  processes.  The  problem  is  to 
mode!  the  essential  characteristics  of  the  Q(z)  profiles  in  a 
way  that  recognizes  operational  observing  and  forecasting 
limitations  yet  takes  maximum  advantage  of  existing  capa¬ 
bilities.  A  prominent  feature  of  the  daytime  aircraft 
soundings  over  inland  areas  is  the  marked  tendency  for 
Q(z)  to  remain  essentially  constant  with  height  within  the 
surface  haze  layer,  and  also  in  the  upper  troposphere 
above  the  primary  haze  layer.  It  should  be  emphasized 
that  the  assumption  of  constant  scattering  ratio  with  height 
does  not  hold  well  for  ground-based  stable  layers  with  lit¬ 
tle  vertical  mixing  such  as  those  associated  with  the  noc¬ 
turnal  formation  of  fog.  However,  for  application  to  prob¬ 
lems  of  contrast  transmittance  in  hazy  atmospheres  in  the 
daytime  following  the  dispersion  of  any  surface  inversion 
existing  at  sunrise,  it  is  reasonable  to  employ  a  model  con¬ 
sisting  of  two  or  more  tropospheric  layers  of  constant  opti¬ 
cal  scattering  ratio.  Thus,  the  forecasting  problem  is 
reduced  to  the  prediction  of  the  altitude  limits  of  the 
atmospheric  layers  and  the  scattering  ratio  within  each 
layer. 

Representation  of  Q(z)  with  a  simple  model  consist¬ 
ing  of  two  tropospheric  layers  was  tested  (Johnson  etal., 
1979)  using  data  obtained  from  experimental  flights  with 
the  instrumented  aircraft  in  western  Europe  in  the  spring 
and  fall  of  1976  and  again  in  the  summer  of  1977.  The 
optical  scattering  ratio  was  assumed  constant  with  altitude 
both  in  the  primary  low-level  haze  layer  and  in  the  rela¬ 
tively  clear  region  aloft  in  the  middle  troposphere.  An 
exponential  decrease  of  scattering  ratio  with  height  was 
assumed  in  the  transition  layer  between  the  iop  of  the 
haze  layer  and  the  base  of  the  upper  troposphere  layer. 
The  depth  of  this  intermediate  transition  layer  was  held 
constant  at  300m.  which  was  the  approximate  average 
transition  zone  depth  for  ihe  measured  profiles.  However, 
the  depth  of  the  low-level  haze  layer  was  made  a  variable 
to  be  specified  or  predicted  in  operational  practice  from 
meteorological  variables.  Attenuation  over  slant  paths 
between  the  surface  and  upper  troposphere  is.  of  course, 
sensitive  to  the  depth  of  the  haze  layer,  /,  as  we’!  ay  to 
the  scattering  ratios  Oi  and  (?„  in  the  primary  (lower!  haze 
la>er  and  in  the  upper  troposphere  layer  respectively. 


Processed  data  from  the  first  three  '''.tiropean  deploy¬ 
ments  for  experimental  flights  reaching  an  altitude  of  at 
least  1!  hm  were  used  for  the  prelim  mar-  tests  of  the 
scattering  ratio  model.  Each  flight  usually  included 
scattering  coefficient  profile  measurements  with  each  of 
four  spectral  filters  having  peaK  wavelengths  of  475,  550, 
660  and  750  nm.  The  four  profiles  were  measured 
sequentially  over  approximately  a  2-hour  period.  Sample 
profile  data  are  shown  in  Fig.  3-i  along  with  the  model 
representations  as  determined  by  the  method  of  least 
squares.  Tne  logarithmic  (base  e)  rms  specification  error 
of  the  individual  profile  values  of  scattering  ratio  meas¬ 
ured  at  30m  intervals  is  15  (or  approximately  15  percent 
rms  error)  for  fiight  C'- 372  and  .1'  lor  flight  C-390.  The 
median  logarithmic  rms  error  of  daytime  profile  represen¬ 
tation  by  the  single  haze  layer  mode!  was  .20  for  the  87 
individual  filter  profiles  (23  flights)  in  the  test  series. 

Limitations  of  the  simple  2-layer  model  representa¬ 
tion  nf  scattering  ratio  profile  result  primarily  from 
occurrences  of  multiple  haze  layer  structure  rather  than 
the  lack  of  scattering  ratio  uniformity  within  a  given  haze 
layer.  For  example,  both  a  dense  haze  layer  near  the 
ground  and  an  intermediate  haze  layer  of  lesser  density 
are  evident  in  the  third  profile  shown  in  Fig.  3-1.  Flight 
C-418  was  an  early  morning  flight  made  over  open  farm¬ 
land  in  northwestern  Germany.  Extension  of  the  model 
to  include  definition  of  the  multi-layered  structure  is 
required  in  such  situations  to  significantly  reduce  the  loga¬ 
rithmic  rms  specification  error  in  this  case,  which  was  near 
0.4  for  the  simple  2-layer  representation  of  the  profile. 
The  computer  code  for  the  calculation  of  contrast 
transmittance  has  provision  for  the  inclusion  of  as  many 
atmospheric  layers  as  warranted  in  complex  situations 
when  detailed  information  about  the  scattering  ratio  profile 
is  available. 


3.2  Specification  of 
Haze  Layer  Depth 

As  summarized  by  Benkley  and  Schulman  (1975), 
several  techniques  have  been  proposed  for  estimating 
boundary  layer  mixing  depths  for  inland  areas.  Various 
models  are  available  for  the  calculation  of  both  the  height 
of  the  nocturnal  boundary  layer  which  is  dependent  pri¬ 
marily  upon  radiational  cooling  processes  and  wind-shear 
turbulence,  and  the  depth  of  the  daytime  mixed  layer 
which  is  dependent  primarily  upon  the  vertical  penetration 
of  free  convection  arising  from  surface  healing.  Particu¬ 
larly  in  the  daytime,  the  ambient  temperature  profile  can 
be  used  to  help  identify  the  expected  mixing  layer  depth 
and  therefore  the  height  of  the  low-level  haze  layer 

Simultaneous  high-resolution  profiles  of  optical 
scattering  ratio  and  temperature  measured  by  the  instru¬ 
mented  aircraft  were  used  to  explore  temperature  lapse 
rate  criteria  for  the  specification  of  haze  layer  depth  (John¬ 
son  etal.,  1979).  The  data  base  consisted  of  the  same 
group  of  23  flights  (87  profiles)  measured  in  western 
Europe  that  were  used  for  preliminary  tests  of  the  scatter¬ 
ing  ratio  model.  Sample  calculations  of  the  temperature 
lapse  rale  profile  and  the  corresponding  scattering  ratio 
profile  are  shown  in  Fig.  3-2.  The  lapse  rate  was  calcu¬ 
lated  for  overlapping  510m  intervals  spaced  30m  apart, 
and  is  plotted  at  the  base  of  the  altitude  interval  over 
which  it  was  calculated.  The  close  relationship  observed  in 
this  case  between  the  base  altitude  of  the  stable  layer  and 
the  depth  of  the  primary  haze  layer  is  typical  of  most 
flights  in  the  test  series. 

Tentative  criteria  were  established  from  the  initial 
tests  to  relate  objectively  the  temperature  lapse  rate  struc¬ 
ture  to  haze  layer  depth.  The  stable  layer  was  defined  as 
the  lowest  altitude  layer  in  which  the  lapse  rate  is  less  than 


Fig  3- !  Scattering  ratio  profiles  as  measured  for  flight  0372  (Soesterberg.  Netherlands!,  flight  0380  (Rodby.  Denmark  I  and  flight  0418 
(Alhom.  Germany)  The  heavy  dashed  line  is  the  objective  model  representation  of  the  profile  The  nephelometer  measurements  were  corrected 
for  the  effects  of  incomplete  purging  of  the  instrument 


Fig.  3-2.  Comparison  of  simultaneous  profiles  of  optical  scattering  ratio,  temperature  and  temperature  lapse  rate 


4.5°C  km~l  (calculated  over  500m  vertical  interval).  The 
altitude  of  the  haze  layer  top  was  assumed  to  be  coin¬ 
cident  with  the  altitude  corresponding  to  the  minimum 
value  of  lapse  rate  within  this  layer,  as  illustrated  in 
Fig.  3-2.  If  no  stable  layer  with  an  identifiable  minimum 
of  less  than  4.5®C  km is  observed,  the  altitude  of  the 
minimum  lapse  rate  observed  below  4  km  in  altitude  was 
selected  as  the  upper  limit  of  the  primary  haze  layer. 
These  criteria  are  straight-forward  except  they  leave 
undefined  the  case  when  a  very  low  inversion  or  stable 
layer  was  present,  and  the  base  of  the  layer  was  below  the 
lowest  altitude  sampled  by  the  aircraft  (usually  100m  to 
200m  above  ground  level). 

The  correspondence  between  the  haze  layer  depth, 
Z,  and  the  estimate  of  Z  given  by  the  temperature  lapse 
rate  criteria  stated  above  is  illustrated  by  one  example  in 
Fig.  3-2.  For  purposes  of  the  initial  tests  (Johnson  etal, 
1979)  the  base  of  the  thermally  stable  layer  was  calculated 
from  the  average  of  the  temperature  profiles  (usually  4) 
measured  by  the  instrumented  aircraft  within  about  a  2 
hour  period.  Similarly,  the  value  of  Z  was  obtained  by  a 
simultaneous  fit  by  the  method  of  least  squares  to  all  of 
the  corresponding  profiles  of  Q(z)  as  derived  from  the 
flight  data  for  the  same  time  period.  For  the  test  sample, 
the  standard  error  of  estimate  of  Z  from  the  lapse  rate  cri¬ 
teria  was  123m.  The  departures  tended  to  be  systematic  in 
that  the  haze  layer  typically  extends  somewhat  above  the 
calculated  base-height  of  the  stable  layer  at  low  z  and 
somewhat  below  the  calculated  height  when  the  haze  layer 
is  relatively  deep.  More  information  on  the  validation  of 
the  profile  modelling  approximations  with  the  experimen¬ 
tal  aircraft  measurements  is  given  in  Johnson  et  at. 
(1979). 

3.3  Summary  of  Proposed  Technique  for  the 
Specification  of  the  Scattering  Coefficient  Profile 
in  the  Troposphere 

The  modelling  techniques  for  estimating  slant-path 
contrast  transmittance  will  accept  as  input  a  much  more 


detailed  representation  of  the  scattering  coefficient  profile 
than  that  given  by  the  simple  2-layer  troposphere  model. 
The  atmosphere  may  be  divided  into  as  many  layers  as 
warranted  by  available  observations  and  forecasts,  with 
each  layer  having  a  different  but  constant  value  of  Q(z). 
However,  a  general  review  of  all  published  high  resolution 
profiles  of  total  volume  scattering  coefficient  and  tempera¬ 
ture  obtained  from  the  many  deployments  of  the  instru¬ 
mented  C-130  aircraft  in  the  United  States  and  Europe 
and  the  results  of  the  limited  quantitative  analysis  based 
on  a  selected  series  of  test  data  suggest  the  following. 

(a)  The  simple  2-layer  troposphere  model  with  a  single 
low-level  haze  layer  of  variable  depth  and  a  rela¬ 
tively  clear  troposphere  aloft  represents  the  vertical 
profile  of  optical  scattering  ratio  with  good  approxi¬ 
mation  in  most  daytime  situations  encountered  by 
the  experimental  aircraft  operating  under  visual 
flight  rules. 

(b)  A  more  complete  haze-layer  model  is  required  for 
reasonable  profile  representation  in  conditions 
identified  by  the  presence  of  a  significant  ground 
based  stable  layer  with  moderate  to  dense  haze  or 
fog  in  addition  to  an  intermediate  layer  of  lower 
haze  density. 

(c)  In  the  absence  of  other  measurements,  the  depth  of 
the  primary  haze  layer  can  be  estimated  effectively 
by  temperature  lapse  rate  criteria  which  identify 
thermally  stable  layers  in  the  lower  troposphere. 
The  relevant  stable  layers  usually  are  of  sufficient 
depth  and  strength  to  be  identified  with  conven¬ 
tional  radiosonde  observations.  However,  standard 
methods  of  processing  the  raw  radiosonde  measure¬ 
ments  result  in  diminished  resolution  of  tempera¬ 
ture  lapse  rate  structure  and,  on  occasion,  seriously 
limit  the  determination  of  mixing  depth. 

(d)  To  a  good  first  approximation,  (he  scattering  mixing 
ratio  can  be  assumed  constant  with  height  in  the  pri¬ 
mary  haze  layer  and  again  in  the  troposphere  above 


Hr  transition  zone.  It  follows,  for  example,  that  the 
average  scattering  ratio  in  the  low-level  haze  layer 
may  be  estimated  effectively  by  an  observation  or 
reliable  forecast  of  scattering  ratio  at  the  base  of  the 
layer. 


4.  SINGLE  SCATTERING  PHASE 
FUNCTION  REPRESENTATION 

As  discussed  in  Section  2,  the  component  of  direc¬ 
tional  path  function  generated  by  the  scattering  of  direct 
solar  radiation  reaching  the  path  is  given  by  the  product  of 
the  phase  function  for  single  scattering,  the  scattering 
coefficient  and  the  solar  irradiance.  The  phase  function 
for  single  particle  scattering  depends  in  a  complex  way  on 
the  size  distribution  and  refractive  index  of  the  aerosol 
particles  and  the  wavelength  of  the  incident  light.  In  gen¬ 
eral,  as  the  particle  size  increases  with  respect  to  the 
wavelength,  the  amount  of  energy  scattered  in  directions 
close  to  that  of  the  incident  radiation  increases  markedly 
while  the  concomitant  backscattering  components 
decrease,  resulting  in  larger  asymmetry  of  the  phase  func¬ 
tion.  Since  the  total  volume  scattering  coefficient  varies 
systematically  as  the  square  or  larger  power  of  the  particle 
radius  depending  upon  the  size  parameter,  we  might 
expect  that  the  total  volume  scattering  coefficient  or 
scattering  ratio,  Q(z),  might  provide,  through  analytic 
representation,  a  good  first  approximation  for  the  single 
scattering  phase  function.  As  discussed  below,  experi¬ 
mental  evidence  shows  this  to  be  true. 

The  approach  used  herein  for  model  development 
was  first  to  represent  the  single  scattering  phase  function 
for  Mie  scattering  by  two  term  Henyey- 

Greenstein  (Irvine,  1968)  functions  as  follows, 

PuiP.g -  c  PHC(p,g{)  +  (1-c)/>„g(0.*2>  ■  (21) 

where 

/W0.*)  -  (1  -  g2)/  I  47t(1  -  2«cos/3  +  g2)}n  I  ,  (22) 
and  the  asymmetry  factor,  g,  is  given  by 

g  -  y  J  P(/3)cos0sin0d0  .  (23) 

The  phase  function  for  single  scattering  has  the  normal¬ 
ized  form. 


The  single  scattering  phase  function  for  combined 
aerosol  and  molecular  scattering  is  bounded  by  the  Ray¬ 
leigh  phase  function  (Q-l,g-Q)  for  a  clear  atmosphere  and 
by  a  phase  function  corresponding  to  dense  clouds  or  fog 
for  very  large  00).  Combined  analysis  of  the  asymmetry 
parameters  (g\,g2 ,c),  as  derived  from  both  a  least-squares 
fit  of  the  Henyey-Greenstein  functions  to  the  average 
phase  functions  measured  by  Barteneva  (1960)  and  to  the 
phase  functions  derived  from  Mie  calculations  using  typi¬ 
cal  aerosol  distributions  in  haze  and  fog,  led  to  the  follow¬ 
ing  consistent  set  of  empirical  equations.  The  equations 
express  the  asymmetry  parameters  as  a  function  of  the 


Mie  scattering  ratio,  (Q-l),  as  follows: 
for  1<Q<1000 

g,(z)  -  .0539  ln[  Q(z)  -  1  I  +  .541  ,  (25) 

g2(z)  -  - .0442  Ini  <?(z)  -  1  ]  -  .471  ,  (26) 

for  OSslOOO 

g,(z)  -  913  ,  (27) 

g2(z)  -  -.776  .  (28) 

for  1  <  @<64.7 

c(z)  -  -  .007  I  ln(p(z)  -  1)  P  (29) 

+  071n|  <?(z)  -  1  ]  +  .81  , 

and  for  Q>64.7 

r(z)  -  0  98  .  (30) 


The  asymmetry  factor,  gM  ,  for  Mie  scattering  is 
given  by 

g/*(z)  -  r(z)sj(z)  +  1 1  -  c(z) )  £j(z) ,  (31) 

and  the  asymmetry  factor  for  combined  aerosol  and 
molecular  scattering  is 

g(z)  -  guU)  [  Q(z)  -  1  MQU) .  (32) 

Similarly,  the  single  phase  function  for  combined 
particle  and  air  molecule  scattering  is  given  by 


P(z,0,4)d n  -  1 


P(fl,z)  - 


VgH-J  Q(z)  -  1 )  PM(0,z) 


In  turn,  the  Henyey-Greenstein  function  parameters, 
;i(z),  g2(z)  and  r(z)  are  approximated  as  individual  func¬ 
tions  of  the  Mie  scattering  ratio,  Q(z)-1,  as  described 
below. 


where  the  theoretical  Rayleigh  phase  function  is 

Pg(fi)  “  3(1  +  cosJ/9)/16ir  .  (34) 


9 


Independent  evidence  of  the  general  applicability  of 
the  model  estimates  of  Pit 3)  with  respect  to  wavelength 
and  altitude  is  given  by  a  comparative  study  by  Johnson 
era/.,  (1979)  of  the  directional  scattering  functions  meas¬ 
ured  by  the  Visibility  Laboratory  airborne  nephelometer 
and  the  Barteneva  (1960)  phase  function  measurements 
In  addition  to  the  measurement  of  total  volume  scattering 
coefficient,  the  integrating  nephelometer  measured, 
separately,  the  directional  scattering  function  at  nominal 
scattering  angles  of  30°  and  150°  in  four  wavelength  bands 
centered  near  475,  550,  660,  and  750  nm.  The  analysis  by 
Johnson  et  a/..  (1979)  of  the  airborne  data  gathered  by  the 
Visibility  Laboratory  in  Europe  and  the  United  States  over 
all  seasons  indicates  that  the  derived  relationship  between 
the  average  phase  functions  and  the  median  scattering 
ratio  for  the  individual  asymmetry  classes  as  measured  by 
Barteneva  (1960)  holds  well  over  all  visible  wavelengths 
and  over  all  altitudes  up  to  the  highest  levels  sampled  by 
the  instrumented  aircraft  (usually  near  6  km). 

In  Fig.  4-1,  the  phase  function  for  single  scattering, 
Pip),  that  was  calculated  using  the  above  system  of  equa¬ 
tions  is  shown  in  comparison  with  sample  calculations 
from  Mie  theory  (Kattawar,  1975)  and  in  comparison  with 
the  Barteneva  measurements  for  a  hazy  atmosphere  with  a 


comparable  scattering  ratio.  The  representation  of  Pip) 
with  the  system  of  two  term  Henyey-Greenstein  functions, 
which  are  in  turn  related  to  Qiz),  are  in  good  agreement 
with  the  Barteneva  function  but  depart  systematical!;,  in 
the  range  of  the  values  extrapolated  by  Barteneva  for 
scattering  angles  outside  the  actual  range  of  the  measure¬ 
ments,  0°  to  16°  and  164°  to  180°.  In  comparison  with 
the  theoretical  phase  functions,  the  model  calculations 
assuming  Qiz)  =  26  (equivalent  surface  visibility  of  about 
10  km)  compare  closely  with  the  phase  function  calcula¬ 
tions  based  upon  the  Haze  L  particle  size  distribution  of 
Deirmendjian  (1969),  and  a  assumed  refractive  index 
m  =  1 .55 — .00/.  The  largest  disparities  are  noted  for 
scattering  angles  near  120°  where  the  model  calculations 
are  higher,  and  near  180°  where  the  model  calculations  are 
lower  than  the  theoretical  calculations. 

A  comparison  of  the  model  and  theoretical  calcula¬ 
tions  of  Pip)  for  a  cloud  case  is  shown  in  Fig.  4-2.  The 
values  of  Pip)  calculated  (Kattawar,  1975)  using  the  indi¬ 
cated  drop  size  distribution  and  refractive  index  are  com¬ 
pared  with  the  Pip)  calculated  from  the  model  assuming 
the  value  of  Q  is  greater  than  1000.  Although  an  exten¬ 
sion  of  model  calculations  to  cloud  situations  yields  rea¬ 
sonable  definition  of  the  general  features  of  the  Pip)  dis- 


Fig.  4-1.  Single  scattering  phase  functions  (a)  calculated  from  Mie  theory 
by  Kattawar  (1475)  using  Diermendjian  (1969)  Haze  L  distribution  with  a 
refractive  index  of  1.55  -  O.Oi;  (b)  observed  by  Beneneva  (I960)  (class  6); 
and  (c)  calculated  from  Eqs.  (25-14)  for  scattering  ratio  of  26 


Fig.  4-2.  Single  scattering  phase  functions  (solid  line)  calculated  from  Mie 
theory  by  Kattawar  (1975)  for  cloud  drop  distribution  n  “  tr6c'' Sr  and 
refractive  index  of  133-  0.0i;  and  (dashed  line)  calculated  from  Eqs 
(25-34)  for  scattering  ratio  of  >  1000. 
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tribution,  important  details  typical  of  phase  functions  for 
fog,  such  as  the  prominent  minimum  near  0-100°  and 
the  secondary  maximum  near  140s,  are  smoothed  out  by 
the  representation  with  Henyey-Greenstein  functions. 

To  the  extent  that  more  accuracy  is  desired  and 
more  complete  information  is  available  to  define  F(z,0), 
the  overall  computer  code  for  calculating  path  radiance 
and  contrast  transmittance  was  made  general  so  that  it  will 
accept  as  input  for  each  atmospheric  layer  any  specified 
P(z,0).  However,  in  the  absence  of  substantive  informa¬ 
tion  other  than  an  estimate  of  Q(z),  calculation  of  P(z,0) 
from  the  system  of  model  equations  is  recommended. 
Additional  evidence  of  the  applicability  of  the  model  for 
estimating  P(z)  from  Q(z)  is  given  by  the  comparisons  of 
measured  and  computed  sky  and  terrain  radiances  that  are 
discussed  in  the  following  section. 


5.  COMPUTER 

application 
OF  MODELLING 
TECHNIQUES 

The  procedures  for  the  calculation  of  sky  and  terrain 
radiances  were  designed  to  take  advantage  of  all  relevant 
input  data.  In  particular,  the  techniques  for  modelling  the 
scattering  coefficient  and  scatter,ng  phase  function  profiles 
and  the  procedures  for  radiative  transfer  calculation  are 
optional  and  may  be  replaced  by  more  detailed  data  or  by 
other  modelling  techniques.  A  summary  of  the  computa¬ 
tional  steps  in  the  assembled  program  and  a  listing  of  the 
input  data  and  the  resultant  products  are  given  in 
Appendix  B. 

As  input  to  the  radiance  and  contrast  transmittance 
computations,  one  may  introduce  as  many  atmospheric 
layers  as  warranted  by  the  completeness  of  the  observa¬ 
tions  and  forecasts.  Minimum  input  data  for  each  layer 
are  the  altitude  limits,  the  average  scattering  ratio  and  the 
average  single  scattering  albedo.  A  listing  of  the  specific 
input  data  for  the  model  calculations  described  below  is 
given  in  Appendix  C. 


5.1  Model  Results 
and  Comparisons 
with  Measured  Data 

The  overall  effectiveness  of  the  modelling  tech¬ 
niques  described  in  Sections  2,  3  and  4  was  tested  through 
direct  comparisons  of  sky  and  terrain  radiances  measured 
bv  the  airborne  scanning  radiometer  and  the  radiance  dis¬ 
tributions  calculated  with  the  model.  A  series  of  compara¬ 
tive  results  are  shown  in  figs.  5-1  -  5-8.  They  illustrate 
some  cf  the  model  performance  characteristics  at  selected 
altitudes  for  several  flight  profiles  covering  a  range  of 
low-level  haze  conditions,  wavelength  intervals  and  solar 
zenith  angles 


Comparative  results  for  flight  C-466,  conducted  over 
a  rural  area  in  northwestern  Germany,  are  shown  in  Figs. 
5-1  and  5-2  for  the  narrow  band  blue  filter  with  a  central 
wavelength  of  475  nm.  This  midday  flight  was  made 
under  conditions  of  isolated  patches  of  middle  clouds  and 
scattered  thin  cirrus.  Moderate  haze  extended  to  1.3  km 
and  the  surface  visibility  was  near  15  km.  The  solar  zen¬ 
ith  angle  averaged  about  44  degrees.  The  9  >  90  s  seg¬ 
ments  of  the  observed  radiance  distribution,  looking 
downward  from  an  altitude  of  200m,  shown  in  Fig.  5-1, 
reflect  significant  fluctuations  that  are  caused  by  changes 
in  the  inherent  reflectance  of  the  underlying  surface, 
which  varied  from  dark  woods  to  open  fields.  As  shown 
in  Fig.  5-2,  the  influence  of  the  surface  reflectance  varia¬ 
tions  on  the  radiance  is  greatly  smoothed  by  the  effects  of 
the  intervening  atmosphere  as  measured  from  an  altitude 
of  6  km  by  the  scanning  radiometer.  The  5-degree  field  of 
view  of  the  radiometer  also  results  in  smoothing.  The 
measured  directional  radiance  values  at  6  km  were  subject 
to  significant  stray  light  error  for  viewing  angles  near  the 
sun  and  are  not  shown  in  the  comparison.  A  prominent 
characteristic  of  the  radiance  distribution  at  6  km  is  the 
bright  horizon  sky  radiance  which  in  this  case  had  compar¬ 
able  brightness  in  the  upsun  (<t>  -  0  °)  and  downsun 
(d>  -  180°)  azimuthal  directions. 


Comparative  results  for  a  Winter  flight  with  low  alti¬ 
tude  sun  conditions  are  shown  in  Figs.  5-3  and  5-4.  Right 
C-401  was  made  in  the  early  afternoon  over  open  farm¬ 
land  and  occasional  woods  in  north  western  France.  A 
very  light  haze  layer  extended  to  an  altitude  of  1400m 
with  a  surface  visibility  of  30  to  40  km.  Scattered  cirrus 
and  low-level  cumulus  clouds  were  present  in  the  area 
during  the  period  of  measurement.  The  results  in  this 
case  are  for  the  broad  band  photopic  filter  with  a  central 
wavelength  of  550  nm.  The  correspondence  in  measured 
and  model  radiance  is  good  at  both  levels  except  for  the 
tendency  to  overestimate  the  radiance  for  zenith  viewing 
angles  near  180s  (directly  downward  below  the  aircraft) 


The  comparative  results  for  a  case  with  thick  haze 
conditions  is  shown  in  Figs.  5-5  -  5-8.  For  flight  C-469, 
the  low-level  haze  layer  extended  to  altitude  1 000m  with 
the  surface  visibility  in  the  range  5  to  8  km.  A  second 
haze  layer  of  somewhat  lower  density  with  the  base  at 
1200m  and  top  1500m  persisted  throughout  the  period  of 
measurement.  It  was  an  afternoon  flight  made  over  the 
central  Netherlands.  The  underlying  terrain  features  were 
predominantly  green  fields  with  occasional  brown  fields, 
woods  and  small  towns.  The  flight  track  was  adjusted 
somewhat  during  the  afternoon  to  remain  in  advance  of  a 
broken  to  overcast  layer  of  cumulus  clouds  located 
immediately  to  the  west  and  north  of  the  track.  Shown 
for  comparison  in  this  case  are  (he  results  for  the  two 
filters  with  central  wavelengths  near  the  upper  and  lower 
limits  of  the  visible  spectrum,  the  blue  filter  (475  nm) 
measure  rtents  are  shown  in  Figs.  5-5  and  5-6,  and  the 
near  infrared  filter  (750  nm)  in  Figs.  5-7  and  5-8 


RADIANCE,  (Wn-'m-W1)  £3  RADIANCE, 
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.  5-3  Comparison  of  measured  and  calculated  sky  and  terrain  radiance  at  altitude  996m  for  night  C-401 
uz,  France).  The  solar  zenith  angle  us  72  deg. 


Fig.  5-4.  Same  as  Fig.  5-3  except  that  the  altitude  is  5290m  and  the  solar  zenith  angle  is  76  deg 
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Fig.  5-5.  Comparison  of  measured  and  calculated  sky  and  terrain  radiance  at  altitude  277m.  Filler  2  (475  nm), 
for  flight  C-469  (Soesterberg,  Netherlands)  The  solar  zenith  angle  is  40  deg 


Fig.  5-6  Same  as  Fig  5-5  except  that  the  altitude  is  5838m  and  the  solar  zenith  angle  is  44  deg 
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Fig.  5-7.  Same  as  Fig  5-5  excepi  it  is  for  Filter  5  (750  nm),  and  the  altitude  is  280m  and  the  solar  zenith  angle 
is  48  deg. 
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Seme  as  Fig.  5-7  except 


6.  SENSITIVITY  OF  SPECTRAL 
CONTRAST  TRANSMITTANCE  TO 
SELECTED  CHANGES  IN  PHYSICAL 
PARAMETERS 

Systematic  model  calculations  were  made  to  illus¬ 
trate  the  effects  of  typical  variations  in  physical  parameters 
upon  the  slant  path  contrast  transmittance  and  the 
corresponding  range  of  object  detection.  For  this  purpose, 
a  two-layer  atmosphere  was  assumed,  consisting  of  a  pri¬ 
mary  haze  layer  of  varying  optical  thickness  and  single 
scattering  albedo,  and  an  overlying  upper  troposphere- 
stratosphere  layer  of  relatively  clear  air  with  a  constant 
optical  scattering  ratio  of  1.3  and  a  single  scattering  albedo 
of  0.97.  The  trial  calculations  were  carried  out  for  an 
assumed  wavelength  of  550  nm.  The  responses  in  con¬ 
trast  transmittance  to  parameter  changes  would  be  qualita¬ 
tively  similar  for  other  visible  wavelengths  but  would  show 
systematic  quantitative  differences  for  the  same  aerosol 
loading  in  the  haze  layer  because  of  the  wavelength 
dependence  of  both  the  total  optical  depth  and  the  single 
scattering  phase  function. 

Results  are  shown  for  a  range  of  haze  layer  depths, 
aerosol  concentrations  and  absorption,  and  for  selected 
values  of  surface  reflectance  and  solar  zenith  angle.  The 
parameters  are  varied  individually  while  holding  each  of 
the  other  variables  constant  and  equal  to  their  values  for 
the  assumed  reference  atmospheres  given  in  Appendix  C, 
which  are  in  general  representative  of  clear  sky  conditions 
observed  in  the  summer  season  in  northern  Germany. 
The  first  two  sets  of  calculations,  Sections  6.1  and  6.2, 
were  based  upon  reference  atmosphere  "A".  The  vertical 
extent  of  the  haze  layer  for  this  reference  atmosphere  was 
assumed  to  be  3  km  except  that  a  few  comparative  calcula¬ 
tions  were  made  also  for  an  assumed  haze  layer  top  at 
500m  altitude  while  holding  the  overall  optical  thickness 
the  same  as  for  the  3  km  case.  The  third  set  of  trial 
results  (Section  6.3)  were  based  upon  reference  atmo¬ 
sphere  ”B"  (Appendix  C),  with  the  haze  layer  extending  to 
1.3  km  and  with  a  single  scattering  albedo  of  .83  and  an 
optical  scattering  ratio  of  16. 

6. 1  Object  at  the  Earth's  Surface 
and  the  Sensor  at  6  Km 

From  Eq.  (4),  the  contrast  transmittance, 
Tc  —  C/C,,  of  a  downward  path  of  sight  to  an  object  at 
the  earth's  surface  is  given  by  the  product  of  the  path 
beam  transmittance,  Tt,  and  ratio  of  the  inherent  back¬ 
ground  radiance  of  the  surrounding  terrain  and  the 
apparent  background  radiance  as  measured  at  the  sensor 
altitude,  bL0lbLr-  It  follows  that  the  contrast  transmit¬ 
tance  decreases  with  an  increase  in  the  optical  thickness  of 
the  viewing  path  in  concert  with  both  a  decrease  in 
Tr(z,9 )  and  an  increase  in  the  path  radiance,  L‘,  gen¬ 
erated  by  the  enhanced  aerosol  particle  scattering.  For  the 
case  with  the  object  at  the  surface  and  the  sensor  at  6  km, 
the  responses  of  Tc  to  changes  in  optical  thickness  and  to 
changes  in  the  cosine  of  the  zenith  viewing  angle  are 
shown  in  Fig.  6-1.  Since  we  have  assumed  that  the  sur¬ 


face  reflectance  obeys  Lambert’s  law,  the  asymmetry  of  Tr 
in  the  upsun  and  downsun  viewing  directions  is  small. 
Except  for  viewing  angles  near  the  horizon,  the  calculated 
Tc  from  the  surface  to  6  km  decreases  more  or  less  in 
direct  proportion  to  the  increase  in  slant  path  distance. 

For  diagnostic  purposes  let  us  define  an  arbitrary 
target  acquisition  range  (TAR)  where  the  contrast 
transmittance  reduces  to  10  percent.  Illustrated  in  Fig.  6-2 
are  the  calculated  changes  in  TAR  that  correspond  to 
selected  changes  in  the  physical  variables.  Again  we 
assume  a  target  located  at  the  earth's  surface  and  a  sensor 
altitude  of  6  km.  The  calculations  were  made  for  an 
azimuth  viewing  angle  of  180°  (downsun).  The  calculated 
TAR  for  the  reference  atmosphere  is  13.5  km.  Comparing 
responses  of  the  TAR  to  the  departures  of  the  individual 
parameters  from  their  values  given  by  the  reference  atmo¬ 
sphere,  we  note  that  a  25  percent  increase  in  T  AR  is  asso¬ 
ciated  with  about  a  25  percent  decrease  in  the  optical 
thickness  of  the  haze  layer.  On  the  other  hand  about  a  50 
percent  increase  in  optical  thickness  is  required  for  a  25 
percent  decrease  in  TAR.  The  dotted  line  shown  in 
Fig.  6-2  for  the  relationship  of  TAR  with  optical  thickness 
depicts  the  calculated  TAR  for  an  haze  layer  top  at  0.5  km 
in  lieu  of  the  3  km  altitude  assumed  for  the  reference 
atmosphere.  Notice  that  the  TAR  differs  only  by  approxi¬ 
mately  15  percent  for  a  haze  layer  with  the  greatly  reduced 
vertical  extent  but  with  increased  extinction  so  as  to  yield 
the  same  optical  thickness. 

Figure  6-2b  shows  that  changes  in  In  TAR  are 
directly  proportional  to  the  changes  in  the  cosine  of  the 
solar  zenith  angle.  For  example,  a  decrease  in  solar  zenith 
angle  of  12  degrees  results  in  about  a  25  percent  increase 
in  TAR  for  the  reference  atmosphere.  Notice  also 
(Fig.  6- 2d)  that  a  25  percent  decrease  in  TAR  is  associated 
with  an  increase  in  haze  layer  absorption  from  the  small 
amount  associated  with  clear  remote  atmospheres  (single 
scattering  albedo  =  0.97)  to  the  large  amount  correspond¬ 
ing  to  urban  atmospheres  (single  scattering 
albedo  =  0.63) 

For  objects  viewed  against  a  terrain  background,  the 
TAR  has  a  significant  direct  dependence  upon  the  surface 
reflectance.  For  example  (Fig  6-2c),  the  calculated  TAR 
decreases  by  25  percent  as  the  surface  reflectance  changes 
from  the  value  associated  with  a  brown  field  (0  10)  to  a 
value  representative  of  a  grass  field  (0.06). 

6.2  Object  at  6  Ft.  and  the 
Sensor  at  Ground  Level 

Diagnostic  model  calculations  of  7j  and  TAR  are 
shown  in  Fig.  6-3  and  6-4  for  an  object  at  6  km  when 
viewed  against  a  sky  background  from  the  earth's  surface. 
The  reference  atmosphere  and  selected  changes  to  the 
reference  atmosphere  are  the  same  as  for  the  sensitivity 
analysis  described  in  the  preceding  paragraphs.  In  contrast 
with  downward  paths  of  sight,  the  calculated  contrast 
transmittance  of  upward  paths  have  large  asymmetry  in 
the  upsun  and  downsun  directions  because  of  the 
enhanced  path  radiance  generated  by  the  increased  scatter- 
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OBJECT  ALTITUDE  -  SURFACE 
SENSOR  ALTITUDE  -  6  km 


80  40  20  I  5  10  8  7  6  7  8  10  15  20  40  80 

SLANT  PATH  DISTANCE  (KM1 


Fit.  6*1.  Diagnostic  calculations  of  contrast  transmittance  as  a  function  of  zenith  viewing  angle  for  a  target 
viewed  against  a  terrain  background.  Values  of  r  gee  iota!  atmosphere  optica)  depth. 


OPTICAL  THICKNESS  COSINE  SOLAR 

SURFACE  TO  SENSOR  ZENITH  ANGLE 

ALTITUDE 


SURFACE 

REFLECTANCE 


SINGLE 

SCATTERING 

ALBEDO 


(a) 


(b) 


(c)  (d) 


Fit  6-2.  Variations  in  laraet  acquisition  ranee  tTr  -  010)  associated  with  departures  of  selected  perimeters 
from  the  aeeumed  reference  atmosphere  vslues.  The  horizontal  dashed  line  denotes  perimeter  values  for  this 
reference  atmosphere  The  horizontal  bars  identify  a  change  of  ±25  percent  from  the  slant  path  distance  where 
Tc  —  0. 10  for  the  reference  atmosphere,  e.g.  13  5  ±3.4  km. 
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SLANT  PATH  DISTANCE  (KM)  WHERE  Tt  -  0  10  |  >  CONTRAST  TRANSMITT ANCF. 


80  40  20  15  10  8  7  6  7  8  10  15  20  40  80 

SLANT  PATH  DISTANCE  (KM) 


Fig.  6-3.  Same  as  Fig.  6-1  except  for  a  target  viewed  against  a  sky  background 


OPTICAL  THICKNESS 

COSINE  SOLAR 

SURFACE 

SINGLF 

SURFACE  TO  SENSOR 

ZENITH  ANGLE 

REFLECTANCE 

SCATTERING 

ALTITUDE 

ALBEDO 

(a) 


(b) 


(c) 


(d) 


Fig.  6-4  St  me  as  Fig.  6-2  except  for  •  urge!  viewed  against  t  sky  background. 
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ing  of  sunlight  at  forward  scattering  angles.  As  shown  in 
Fig.  6-3,  the  calculated  Tc  for  <t>  -  O'is  less  than  10  per¬ 
cent  for  all  zenith  angles  when  the  total  optical  depth  is 
0.5  or  larger.  In  the  downsun  direction  (<6-180"),  Tc 
has  only  a  modest  variation  with  zenith  viewing  angle  for 
angles  less  than  about  70°. 

A  strong  dependence  of  target  acquisition  range  in 
the  downsun  direction  upon  the  optical  thickness  of  the 
haze  layer  is  evident  in  the  results  shown  in  Fig.  6-4a.  A 
change  in  thickness  of  only  10  percent,  holding  other 
parameters  of  the  reference  atmosphere  constant,  results 
in  about  a  25  percent  change  in  TAR.  A  marked  sensi¬ 
tivity  of  the  TAR  to  changes  in  solar  zenith  angle  is  also 
found,  (Fig.  6-4d).  The  observed  decrease  in  TAR  with 
an  increase  in  solar  altitude  stems  from  the  associated 
increase  in  the  path  radiance.  In  comparison  with  the 
results  shown  in  Fig.  6-2,  the  TAR  for  objects  viewed 
against  a  sky  background  show  significantly  less  response 
to  changes  in  the  single  scattering  albedo  (Fig.  6-4d)  of 
the  haze  layer  and  changes  in  surface  reflectance  (Fig.  6- 
4c). 

6.3  Object  at  the  Surface  and 
the  Sensor  at  20  Km 

The  third  and  last  set  of  trial  calculations  differs 
from  the  previous  two  in  that  a  different  reference  atmo¬ 
sphere  (REF-B  in  Appendix  C)  was  assumed.  The  haze 
layer  is  more  shallow,  extending  to  1.3  km,  and  the 
scattering  ratio  and  single  scattering  albedo  in  the  haze 
layer  are  16  and  0.83  respectively.  With  the  assumed  sen¬ 
sor  altitude  of  20  km,  the  viewing  path  has  a  long  segment 
through  the  relatively  clear  region  in  the  lower  strato¬ 


sphere  and  upper  troposphere.  Yet  another  departure 
from  the  previous  calculations  is  that  the  target  at  ground 
level  is  viewed  in  the  upsun  rather  than  the  downsun 
direction. 

The  calculated  slant  range  corresponding  to  5  per¬ 
cent  apparent  contrast  (contrast  transmittance  of  10  per¬ 
cent)  for  the  reference  atmosphere  and  for  specific 
changes  in  the  individual  parameters  of  the  atmosphere 
are  shown  in  Fig.  6-5.  We  note,  for  example,  that  each  of 
the  following  parametric  changes  in  the  reference  atmo¬ 
sphere  result  in  a  decrease  in  the  slant  path  distance 
corresponding  to  5  percent  contrast  transmittance  in  the 
range  of  19  to  25  percent: 

a.  50  percent  increase  in  haze  layer  depth 

b.  Factor  of  5  increase  in  the  aerosol  scattering  ratio  in 
the  upper  troposphere 

c.  Increase  in  haze  layer  absorption  corresponding  to  a 
change  in  single  scattering  albedo  from  .83  to  .63 

d.  Increase  in  solai  zenith  angle  from  50  to  60  deg 

It  should  be  emphasized  that  the  results  described  in  this 
section  are  specific  only  for  the  conditions  imposed  on  the 
trial  calculations.  However,  these  examples  provide  some 
insight  into  the  changes  in  image  transmittance  charac¬ 
teristics  that  are  associated  with  typical  changes  in  the 
environment. 


7.  SUMMARY  COMMENTS 

The  techniques  developed  during  the  course  of  the 
experimental  optical  measurement  and  analysis  program 


OBJECT  AT  SURFACE  INHERENT  CONTRAST  -  SON 
OBSERVATION  LEVEL  -  23  KM  A7.IMI  TH  VIEW  ANGLE  -  <1* 
WAVELENGTH  -  550  om  (PHOTOPK  > 
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Fig.  6-5.  Comparison  of  the  slam  range  corresponding  lo  apperem  contras! 
of  5  percem  for  reference  atmosphere  (Ref  B,  Appendix  C)  and  selected 
changes  in  ihe  reference  atmosphere 


can  provide  relatively  fast  and  consistent  estimates  of  opti¬ 
cal  atmospheric  properties  over  any  slant  path,  as  a  func¬ 
tion  of  wavelength  in  the  visible  spectrum.  To  the  extent 
that  a  climatological  data  base  of  conventional  meteorolog¬ 
ical  observations  exists,  the  techniques  in  their  present 
form  can  be  used  to  estimate  readily  the  frequency  distri¬ 
butions  of  such  quantities  as  spectral  contrast  transmit¬ 
tance,  as  a  function  of  location  and  season. 

Refinements  in  the  modelling  approximations  are 
being  made  as  validation  tests  and  experiments  continue. 
Results  to  date  indicate  that  further  parameterization  of 
selected  components  of  the  model  can  be  made  which  will 
improve  computational  efficiency  and  at  the  same  time 
retain  the  capability  to  take  advantage  of  all  relevant 
observing  and  forecasting  information.  A  realistic  goal  is 
to  employ  the  computer  model  as  part  of  a  data  acquisition 
and  microprocessing  system  for  real  time  estimates  of 
image  transmission  properties.  Additional  experiments  are 
planned  to  examine  further  the  tradeoff's  between  the 
type,  accuracy,  frequency  and  density  of 
optical/meteorological  observations  and  the  reliability  of 
the  estimates  of  spectral  contrast  transmittance. 
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APPENDIX  A 


and 


DERIVATION  OF  THE  COMPONENT  OF 
PATH  RADIANCE  DUE  TO  DIRECTIONAL 
SCATTERING  OF  DIFFUSE  RADIANCE 
USING  THE  DELTA-EDDINGTON 
APPROXIMATION 
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-l 


(l-f')(l+3g‘nn')dfi‘  -  2(1-/')  , 


Eq.  (A.2)  reduces  to 


(A.4) 


The  following  derivation  of  Eq.  (13)  in  Section  2.1 
was  provided  by  Mr.  Eric  P.  Shettle  of  the  Optical  Physics 
Laboratory,  Air  Force  Geophysics  Laboratory.  The  equa¬ 
tions  referenced  in  the  derivation  are  designated  by 
number  as  given  in  Section  2  of  this  paper. 

Substituting  Eq.  (12)  into  the  last  term  of  Eq.  (6) 
we  have, 


s (z)  f  J*  LDh)PDh,0)eifi'tl6' -  s(2)L0U) .  (A. 5) 

0  -1 


Using  Eq.  (7)  in  the  last  term  of  Eq.  (A.l)  and  integrating 
over  it'  we  have 
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(A. 6) 
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n'l  2f'S(ji-ii')  +  (!-/') (1+3* 'mm  Idu' . 


where  n'  -  cosfl',  and  diY  -  dn’d6'.  Using  Eq.  (7)  in  the 
first  term  on  the  right  hand  side  of  Eq.  (A.l)  and  integrat¬ 
ing  over  <t'  we  have. 
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Where  cos/3  -  mm'  +  sintoinfl'cos(</>-<4') 


and  integrating  over  m’  we  have 
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(A. 7) 


Substituting  Eqs.  (AS)  and  (A.7)  into  Eq.  (A.7)  and  from 
Eq.  (II),  I  /’  +  (1-/’)*’  -  g  ),  the  expression  for  the 
diffuse  component  of  the  path  function  is  given  by 


and 
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-  s(z)  I  Lp(z)  +  gLp  (z) cos»  I  .  (A. 8) 

Also: 

(Eq.  13,  Section  2.1) 
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APPENDIX  B 


SUMMARY  OF  INPUT  DATA, 
COMPUTATION  STEPS  AND 
RESULTANT  DATA 


B.l  General  Input  Data 

B.  11  Number  of  atmospheric  layers 

B.1.2  Representative  wavelength  (pm) 

B.l. 3  Extraterrestrial  solar  (lunar)  irradiance 

( w/m2pm ) 

B.l. 4  Rayleigh  optical  depth  for  top  layer  for 

A  —  ,55pm  (unless  otherwise  specified,  set  equal 
to  .0258  for  layer  base  altitude  of  10  km). 

B.  1.5  Solar  zenith  angles  (deg.) 

B.l. 6  Zenith  viewing  angles  (deg.) 

B.1.7  Scattering  angles  between  viewing  path  and  sun 
(deg.,  max.  10). 

B.  1.8  Average  reflectance  of  underlying  surface  and 
reflectance  of  background  in  the  vicinity  of  target 
if  different  than  the  average  reflectance. 


B.2  Input  Data  Each  Layer 

B.2.1  Base  altitude  (km) 

B.2.2  Scattering  mixing  ratio 

B.2. 3  Single  scattering  albedo 

B.3  Computation  Steps 

B.3.1  Compute  Rayleigh  atmosphere  optical  thickness 
for  each  layer 

B.3. 2  Compute  total  optical  thickness  for  each  layer. 

B.3. 3  Compute  optical  depth  for  the  base  of  each  layer. 

B.3. 4  Compute  single  scattering  phase  function  and 

asymmetry  factors  for  each  layer.  (These  param¬ 
eters  may  be  specified  if  desired  in  lieu  of  the  use 
of  model  algorithms). 

B.3. 5  Transform  optical  depth,  single  scattering  albedo, 

and  asymmetry  factor  in  accordance  with  delta- 
Eddington  approximation. 

B.3.6  Compute  path  function  distribution  for  the  base 
of  each  sub-layer  (each  layer  is  divided  equally 
into  2  sub-layers). 

B.3. 7  Compute  path  radiance  distribution  for  each  sub¬ 

layer. 


B.3. 8  Compute  path  transmittance  distribution  for  each 

sub-layer. 

B.3.9  Compute  sky/terrain  radiance  distribution  at 

sub-layer  boundaries. 

B.3. 10  Compute  slant  path  contrast  transmittance  of 
selected  paths. 


B.4  Output  Tables  and  Data 

B.4.1  Single  scattering  phase  function,  />(z,/ 3),  as  calcu¬ 
lated  or  specified  for  each  layer. 

B.4.2  Path  radiance,  Z.'(z,0,0),  for  each  sub-layer  and 
each  solar  zenith  angle,  0,. 

B.4. 3  Path  transmittance  T,(z,0)  for  each  sub-layer. 

B.4.4  Sky  radiance,  Z.,(z,0,/3),  for  O°$0<85°  for 

selected  levels  (max  5)  for  each  0,. 

B.4.5  Sky  radiance,  L,(z,0,/3)  for  95°<0<18O°  for 
selected  levels  (max  5)  for  each  0,. 

B.4. 6  Sky  plus  terrain  radiance  L,(z,0,0)  for 

95°^0<  180°  for  selected  levels  (max  5)  for  each 

0,. 

B.4.7  Contrast  transmittance,  Tr(z,0,/3)  between  object 
and  sensor  altitudes  (max.  5  altitude  intervals) 
for  0°< 0^85 “and  for  each  0,. 

B.4.8  Contrast  transmittance  Tf  (z,0,/3)  between  object 
and  sensor  altitudes  (max.  5)  for  95°<0<18O° 
and  for  each  0S. 


B.5  Additional  Output  From 
The  Supplemental  Eddington 
Computer  Program 

The  standard  Eddington  computer  program  is 
employed  as  an  integrated  part  of  the  model  calculations 
of  contrast  transmittance.  The  following  output  from  the 
program  as  modified  by  the  delta- Eddington  transforma¬ 
tion  of  input  parameters  is  available  in  tabular  form  for 
the  base  of  each  sub-layer  and  for  each  solar  zenith  angle 

optical  depths,  r  and  r' 

components  of  Eddington  diffuse  radiance.  Lt,  (z> 

and  Lp  i:) 

total  scalar  irradiance 

downward  diffuse  irradiance 

upward  diffuse  irradiance 

solar  irradiance 

total  downward  irradiance 


APPENDIX  C 


SPECIFIC  INPUT  DATA  FOR 
MODEL  CALCULATIONS 


Flight  Number 

Wavelength 

Number  of  Layers 

Tropopause  Altitude 

Upper  Altitude  Limit 

C-401 

550  nm 

3 

10  km 

C-466 

475  nm 

3 

10  km 

C-469 

475  nm 

5 

10  km 

C-469 

750  nm 

5 

10  km 

REF-A 

550  nm 

3 

10  km 

REF-B 

550  nm 

3 

10  km 

3rd  Haze  Layer 

1.8  km 

1.8  km 

2nd  Haze  Layer 

1.5  km 

1.5  km 

1st  Haze  Layer 

1  km 

1.3  km 

1.2  km 

1.2  km 

3  km 

1.3  km 

Optical  Scat.  Ratio 

Stratosphere 

1.2 

1.2 

1.2 

1.8 

1.3 

1.3 

Upper  Trop. 

1.2 

1.2 

1.2 

1.8 

1.3 

1.3 

3rd  Haze  Layer 

7.3 

20.8 

2nd  Haze  Layer 

1.3 

7.0 

1st  Haze  Layer 

5.0 

12.0 

28.0 

54.0 

13.0 

16.0 

Single  Scat.  Albedo 

Stratosphere 

.99 

.99 

.99 

.99 

.97 

.97 

Upper  Trop. 

.97 

.97 

.97 

.97 

.97 

.97 

3rd  Haze  Layer 

.83 

.83 

2nd  Haze  Layer 

.97 

.97 

1st  Haze  Layer 

.97 

.83 

.83 

.83 

.97 

.83 

Surface  Reflectance 

1 

1  .09 

.07 

.07 

.33 

.10 

.07 
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GLOSSARY  AND  NOTATION 


The  notation  used  in  reports  and  journal  articles 
produced  by  the  Visibility  Laboratory  staff  follows,  in  gen¬ 
eral,  the  rules  set  forth  in  pages  499  and  500,  Duntley 
et  al.  (1957).  These  rules  are: 

Each  optical  property  is  indicated  by  a  basic  (parent)  sym¬ 
bol. 

A  presubscript  may  be  used  with  the  parent  symbol  as  an 
identifier,  e.g.,  b  indicates  background  while  t  denotes  an 
object,  i.e.  target. 

A  postsubscript  may  be  used  to  indicate  the  length  of  a 
path  of  sight,  e.g.,  r  denotes  an  apparent  property  as 
measured  at  the  end  of  a  path  of  sight  of  length  r,  while  o 
denotes  an  inherent  property  based  on  the  hypothetical 
concept  of  a  photometer  located  at  zero  distance  from  an 
object,  i.e.  target. 

A  postsuperscript*  or  postsubscript*,  is  employed  as  a 
mnemonic  symbol  signifying  that  the  radiometric  quantity 
has  been  generated  by  the  scattering  of  ambient  light 
reaching  the  path  from  all  directions. 

The  parenthetical  attachments  to  the  parent  symbol  denote 
altitude  and  direction.  The  letter  z  indicates  altitude  in 
general;  z,  is  used  to  specify  the  altitude  of  a  target.  The 
direction  of  a  path  of  sight  is  specified  by  the  zenith  angle 
0  and  the  azimuth  4>.  In  the  case  of  irradiances,  the 
downwelling  irradiance  is  designated  by  d,  the  upwelling 
by  u. 

The  radiometric  symbols  used  herein  now  correspond  to 
the  OSA  recommendations  in  Section  1  of  the  Handbook 
of  Optics,  McGraw-Hill  Book  Co.,  Driscoll  and  Vaughn 
(1978).  Prior  to  June  1980,  the  symbol  used  for  radiance 
l  was  #f,  for  irradiance  e  was  h,  and  for  attenuation  length 
was  /. 


Albedo  al  altitude 
/4<z>-  E(z.u)/E(z,d). 


Limb  darkening  factor  relating  the  average  tun  radi¬ 
ance  and  the  center  sun  radiance  £><x )  -  $L0I SLQ 


Spectral  irradiance  (formerly  symbol  H)  defined  as 
Ek  -  j  icosMfl 


broad  band  sensor  irradiance,  defined  as 


Irradiance  produced  by  downwelling  flux  as  deter¬ 
mined  on  a  horizontal  flat  plate  at  altitude  z  (form¬ 
erly  A(x.d)).  In  this  report  d  is  used  in  place  of  the 
minus  sign  in  the  notation  (//(;,-))  which  appears 
in  Duntley  (1969).  This  property  may  be  defined  by 

the  equation  Eiz,d)  -  J  L(z,8\4')cas$'dt\' 


Irradiance  produced  by  the  upwelling  flux  as  deter¬ 
mined  on  a  horizontal  flat  plane  at  altitude  z  (form¬ 
erly-  Here  u  is  substituted  for  the  plus  Sign 

formerly  used  in  the  notation  (//(z,+)>. 


Asymmetry  factor  J*  P(z,fikosfid(cosfi) 


Scale  height  at  altitude  z,  the  height  of  a  homogene¬ 
ous  atmosphere  having  the  density  of  the  layer  at 
altitude  z. 


Spectral  radiance  (former  symbol  N) 


Broad  band  sensor  radiance  is  defined  as 


L  -  I  L,Sk7\  dk/ok 


Inherent  radiance  based  on  the  hypothetical  concept 
of  a  photometf  -  located  at  zero  distance  from  an 
object  at  altitude  z  in  the  direction  specified  by  zen¬ 
ith  angle  %  and  azimuth  4. 


Apparent  radiance  as  determined  at  altitude  z,  from 
the  end  of  a  path  of  sight  of  length  t  at  zenith  angle 
6  and  azimuth  4-  This  property  may  be  defined  by 
“  L0(zr$.4irf(z,$l  +  l'(z,9,4) 


C,  (!,#.♦) 


Absorption  coefficient 


Inherent  contract  determined  for  e  ptth  of  etpht  of 
tern  lenfth  at  the  altitude  of  the  object  <(  in  the 
direction  of  aenlth  antte  *  and  aaimuih  * 


Apparent  contract  ae  determined  at  altitude  i  from 
the  end  of  path  of  ai|ht  of  length  r  in  the  direction 
of  zenith  anffie  •  and  aaimuth  * 

Radiance  distribution  function 

0<z>  -  f tif.r.dtaacren/tir) 


Stock  body  radiance  at  wavelength  A  and  tempera¬ 
ture  T 


Equilibrium  radiance  at  altitude  z  with  the  direction 
of  the  path  of  sight  specified  by  zenith  angle  f  and 
azimuth  4 ■  lids  property  is  a  point  function  of  posi¬ 
tion  and  direction 


Nth  function  at  altitude  z  with  the  direction  of  the 
path  of  sight  specified  by  zenith  angle  •  and  attmuth 
4.  This  property  is  defined  by  the  aquation 

t.'(z.d.*)-  j  oii.z'Uii.r.dien' 
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f.’  (Z.S,0» 


Wlvm*  ptih  rwjian<e  *5  determined  et  altitude  r  it  the  end 
of  a  path  of  »*ht  of  length  r  in  the  direction  of  icn- 
ith  anile  9  irvd  azimuth 


W/snrr  sky  radiance  «  altitude  z .  zenith  angle  9  and 
azimuth  *  Aleo  the  peth  redience  for  the  peth  of 
eilhi  of  length  °»  from  out  of  the  etmosphere  to  elti- 
tude  x. 


Symbol 

r.U.1) 

w(z) 


W/srm1 


kg/m 1 


Apparent  radiance  of  the  center  of  the  solar  disk  a* 
determined  from  the  end  of  a  path  of  si|ht  of  len|th 
«  from  out  of  the  atmosphere  to  altitude  z  at  the 
zenith  angle  of  the  sun  0j. 


Absolute  air  masa 


it  mile  9  (z.0)  -  J*  p(z)dr. 


a(z) 


Units 


none  Radiance  transmitaftce  for  the  path  of  sifht  at  zenith 

angle  9  from  out  of  the  atmosphere  to  the  altitude  z 


Single  scattering  albedo  or(z)-z(z )/<*(/) 


Altitude,  usually  used  as  above  ground  level 


f*  Altitude  of  any  applicable  target. 

m“l  Volume  attenuation  coefficient  as  determined  at  alti¬ 

tude  z.  a(z)— «(z)+*<*) 

deg  Symbol  for  scattering  angle  of  flux  from  a  light 

source.  It  is  equal  to  the  angle  between  the  line 
from  the  source  to  any  unit  scattering  volume  and 
the  path  of  a  ray  scattered  off  this  direct  line 


none  Relative  optical  air  mass 


A 


m(Pj) 


*(z) 


none  Relative  optical  airmasa  at  zenith  angle  of  the  sun. 
shorthand  for 

<(z) 

none  Refractive  index 


PUS) 


QU) 


Normalized  phase  function  for  single  scattering. 
PUS)  -  <r(z,0)/z<z) 


Optical  scattering  mixing  ratio  at  altitude  z.  This 
quantity  is  defined  as  the  ratio  of  the  total  volume 
scattering  coefficient  at  altitude  z.  to  the  molecular 
(or  Rayleigh)  volume  scattering  coefficient  at  the 
same  altitude  z.  (?(z)-s(z)/^s(z) 


.€<*> 


none  Symbol  to  indicate  incremental  quantity  and  used 

with  r  and  z  to  indicate  small,  discrete  increments  in 
path  length  r  and  altitude  z. 

W/m ^  Srelfr  irradiance  This  may  be  defined  as  the  radiant 

flux  arriving  at  a  point,  from  all  directions  about  that 
point,  at  altitude  z  CTyler  and  i»reisendorfer  (1962)) 

,<z)-  f  £<*,». *><*« ' 


Diffuse  scslsr  irradiance  ,«(*). 


Wfrn^  sun  icMltr  irradiance  al  eliilude  J 

J.tz  )-,«<«) 


l(z) 


„  Path  length,  foe  paths  of  sight  st  zenith  angles  0  to  { 

70  degrees.  r-ttclbt 

„-l  Toul  volume  scstler ng  coefficient  es  determined  st 

altitude  I  This  property  may  he  deffned  by  the 

equations  t(z)  -  f  <e(zd»dfl  -  +  w’*** 


In  the  abeence  of  atmoephehc  ihsorption.  the  total 
volume  scattering  coefficient  is  numerially  equal  to 
the  attenuation  coefficient 


Volume  saiieftng  coefficient  foe  Mie  It  serosol. 
sauenng  si  siniude  z 

p(z) 


m  Radius  of  the  earth. 

deg  Symbol  for  zenith  angle  Thu  symbol  i,  usually  used 

as  one  of  two  coordinates  to  specify  the  direction  of 
a  path  of  sight. 

deg  Symbol  for  zenith  angle  usually  used  as  one  of  iwo 

coordinates  to  specify  the  direction  of  a  discrete  pot- 
lion  of  the  sky. 

deg  Zenith  ingle  of  llte  sun 

nm  Symbol  for  wavelength 


lm/m 3  Density  it  altitude  z 


M~  * 

Volume  sailer  ng  coefficient  for  Raytagh  It  molec¬ 
ular  sal  ter  ng  st  altitude  r 

«lz,») 

r<z> 

•jr 

r,(z.») 

none 

nut  inn  iransntnwta  as  determined  at  altitude  > 
for  a  peth  of  atffffit  of  len*h  r  at  zenith  angle  * 
(formetly  referred  to  ■  "beent-  irammlttanee)  Till* 
property  ta  ntdepervSant  of  azimuth  in  at.vwaidtaroa 
having  horizontal  uniformity  It  it  alweys  tho  aente 
for  the  dangnarad  path  of  ugtit  or  ita  toetproal 

r(z) 

♦ 

Contrast  tranamttance 

r/z.a.o)  -  c„i z,.».*>/c,(..#.*> 

Symbol  for  volume  scntiering  function  Peranthett- 
ctl  symbols  «re  r  to  designate  altitude  end  fl  to 
deMnate  the  scatter f*  angle  from  t  source 


Optical  dspth  r(z)  -  J  n(r>* 
z 


Symbol  for  azimuth  Tht  azimuih  ti  tha  sngla  in  tha 
horizontal  (dm  of  the  Obsarver  between  .  feted 
point  and  the  peth  of  U0u  Tbs  ffsed  pram  mey  h. 
for  eumtM.  tree  North,  the  beano*  of  the  sun.  or 
tht  Peering  of  the  moon  TV*  symbol  *  vmi*r 
used  ■  out  of  taro  coordffiaaet  »  tpadfy  the  dime- 
lion  of  ■  peth  of  u#n 
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limits  SynM 


4‘  dec  This  symbol  for  uimuth  is  usually  used  as  one  or  * 

two  coordinates  to  specify  the  direction  of  a  discrete 
portion  of  the  sky 

n 

4  def  Angular  solar  radius  at  true  earth-to-sun  distance 


limits 

deg 

Angular  solar  radius  at  mean  solar  distance 

sr 

Sym>  *- 

solid  angle 

For  a  hemisphere 

7or  a  sphere: 

0—4 it  steradians 

I 


! 

I 

i 


F 
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